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ABSTRACT 

A general formalism to include experimental reaction cross sections into calculations 
of stellar rates is presented. It also allows to assess the maximally possible reduction 
of uncertainties in the stellar rates by experiments. As an e xample for the a pplica tion 



of the procedure, stellar neutron capture reactivities from iKADoNiS vO.a (|2009l ) are 
revised and the remaining uncertainties shown. Many of the uncertainties in the stellar 
rates are larger than those obtained experimentally. This has important consequences 
for s-process models and the interpretation of meteoritic data because it allows the rates 
of some reactions to vary within a larger range than previously assumed. 

Subject headings: Nuclear reactions, nucleosynthesis, abundances 



1. Introduction 

An increasing number of reaction cross section measurements are performed with the aim 
to improve reaction rates for nucleosynthesis studies. The questions of how to assess the impact 
of the experimentally obtained cross sections on astrophysical reaction rates and how to convert 
laboratory cross sections to stellar rates inevitably arise. Closely related is the question concerning 
the estimate of remaining uncertainties in the stellar rates. These questions are addressed in the 
following. 

The stellar enhancement factor (SEF) - the (theoretically predicted) ratio of stellar and lab- 
or atory rate - was used in the past to derive the stellar rate from a measurement. It was shown 



m 



Rauscher et al.l (|201ll ) that the SEF is not adequate for this purpose and the ground state 



contribution Xq was introduced to replace the SEF. 

Here, a further generalization of this concept is presented and a formalism for including data 
and their uncertainties into an improved stellar rate with revised uncertainties is laid out. The 
generalized relative level contribution to the stellar rate is introduced in § [2j The proper inclusion 
of data into stellar rates is discussed in § [3l The following § |1] then explains how the attached 
error changes when using experimental information to improve a stellar rate. As an important 
application, neutron capture rates for the s-process are revised in § [5j They are commonly believed 
to be strongly constrained by precise data but it will be shown that the remaining uncertainties 
are not dominated by the experimental errors in many cases. 
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The relative contribution X to the stellar rate 



The stellar rate r* is a weighted sum over reactivities Ri of thermally excited states, which are 
each bombarded by Maxwell-Boltzmann energy-distributed projectiles, 



(1) 



with the number den sities Unrn, ^ ta.r of projectile and target nucleus, respectively. The statistical 
weights are given by (jFowleiill974l ) 



Wi 



(2J, + l)exp(-£;f"7(fcT)) 
E^{(2J^ + l)exp(-ii;-7(A;r))} 



(2J, + l)exp(-£;f"7(A;T)) 
G 



(2) 



where Ji and Ef^'^ are spin and energy of the ith target level and G is the nuclear partition 
function, summing over all excited target levels. Each reactivity is derived by folding the reaction 
cross section with the energy distribution $ of the projectiles. 



Ri 



ai{Ei)^{Ei,T)dEi 



(3) 



where in the case of interacting nuclei at usual stellar plasma conditions ^ is the Maxwell-Boltzmann 
distribution. 



The quantity 



-/{kT) 



fai{E)<^{E,T)dE 
J a''^{E)<i>{E,T)dE 



(4) 



describes the contribu t ion of transition s commencing on level i of a target nucleus to the stellar rate 
(jRauscher et al.l l201ll : iRauscherl |2012| ) . The latter also includes additional transitions originating 
from further thermally populated excited states. The reaction cross section cjj is defined as 



and the effective cross section as (|Fowlerlll97i iHolmes et allllQTd ) 

EE' 



.eflf 



(5) 



(6) 



where o"*"^-^ is the partial cross section from target level i to final level j (cross sections at zero or 
negative energies are zero). 

The relative contribution Xi can be used to determine the impact of laboratory measurements 
on stellar rates. Laboratory cross sections only consider reactions on nuclei in a single level i = i^^^. 
Usually, this is the ground state (g.s.) of the nucleus, which implies i^'"^^ = (when counting of the 
levels starts with at the g.s.) and thus a^^^ = Y^jO'^^-'- A few nuclides naturally occur in an 
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isomeric state above the g.s. In this case, the laboratory cross section is derived using the isomeric 

target state iiso and thus 

^lab — Zjgo. Among s-process nuchdcs there is the isomer ''^^^^Ta which is 
a level at Ef^ = 77 keV excitation energy and with spin Jjj^^ = 9. 

Often used in experimental s-process studies, the MACS (Maxwellian Averaged Cross Section) 
(a) is nothing else than the laboratory reactivity i^jiab divided by the most probable projectile 
energy, 

W)iT) = ^ , (7) 

where E = kT is the most probable energy in the Maxwcll-Boltzmann energy distribution. For 
use in nucleosynthesis calculations, a MACS always has to be converted to a rate. Obviously, this 
cannot be a stellar rate r* when < 1. Therefore a rate directly computed from a (measured) 
MACS is only useful when X ^ 1. Otherwise, the rate has to be calculated from theory, perhaps 
using some information from the MACS, as explained in the following. 



3. Renormalization of theory to experiment 
3.1. General 

A theoretical stellar rate r^j^ cannot be directly renormalized to an experimental value because 

■lab 

the measurement only yields a rate on a single level r'^p or MACS {o')^^^. Therefore, only the 
calculated rate rf^^ on level i^'^^ or MACS (<T)th can be compared to the measured one, giving a 
ratio 

,-lab 



(a) 

I ' exp \ / 



V - —Kb - 7-T , (oj 



respectively. But what is actually needed for a renormalization of the stellar rate is the ratio 

v*'='-^ , (9) 

where r*g^ is the new stellar rate including the experimental information. There are two possible 
approaches for implementing the experimental results into a stellar quantity: 

1. Assume that the transitions from the excited states of the target have the same deviation as 
the one found for the transitions from the laboratory level i^'^^; then v*' = v' . 

2. Realize that the measurement of the laboratory level transitions did not provide any infor- 
mation about the transitions originating from other excited states of the target. Therefore 
only the laboratory level contribution to the stellar rate should be renormalized and the other 
contributions left unchanged. This leads to a renormalization factor 

i;*' = f;%iab + (l-X^iab) = l + X^iab(t;'-l) , (10) 

with Xjiab being the laboratory contribution as before. 
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Note that the renormahzation factors are temperature-dependent in both approaches. 

Assumption 1 can be used when the reason for the deviation of theory from experiment is 
known and also that it apphes (with the same magnitude) to the cross sections of the excited levels. 
This cannot be inferred generally but would have to be investigated for each nucleus and reaction 
separately. Often the errors in predicting the transitions from the g.s. and from excited levels may 
be correlated. Nevertheless, it is not easily justified to claim that they are the same. Assumption 
2 presents a more conservative view to the other extreme, implying that the measurement did not 
constrain the cross sections of the other levels in any way and that therefore these should not be 
included in a renormahzation. Only the transitions on the laboratory level are renormalized in that 
picture. 

To be sure, using assumption 2 combined with the error estimate described in §|1] below is the 
safest way to perform a renormahzation. 



3.2. Selection of renormalization temperature 

It is also important to consider the applicability of the theoretical model to calculate the cross 
sections or reaction rates. For compound nuclei with low level density at the compound formation 
energy (for exampl e, nuclei at magic neutron numbers), t he statistical Hauser-Feshbach model may 



not be applicable (jRauscher. Thielemann. &: Kratdll997l ). If this is the case, then also the energy 



dependence will not be correct. The latter is also true when the model is applicable at the higher 
energy part of the considered (n,7) energies but not at lower energies. 

It would be best to use renormalized Hauser-Feshbach values only in the region where the 
model is applicable and employ some other temperature/energy dependence outside that region. 
Unfortunately this is only possible when there are experimental data at the low energies, which 
often is not the case. If experimental data were available, a renormalization would not be required, 
anyway. 

When renormalizing, it has to be taken care whether the renormalization value really is at 
an energy or stellar temperature in the applicability range of the model. If not, a value at higher 
energy or temperature has to be used. Again, this is only possible when experimental data are 
available not only at one energy or temperature. 



4. Determining uncertainty factors and error bars on stellar rates 

Let us assume the theoretical stellar rate r* (or stellar reactivity) has an uncertainty factor 
U* = C/th > 1- For example, a 50% uncertainty translates into C/th = 1-5 and the true rate r^J.^p is 
expected to be in the range r*/1.5 < ^^J.^g < 1.5r*. Let us further assume that an experiment has 
measured the MACS or g.s. rate {i^^^ = 0) with an uncertainty factor f/cxp ^ 1 which is smaller 
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than f7th- When the theoretical MACS is replaced by its experimental value, this will r esult in a 



new stellar rate as described in § [3l The new uncertainty in this stellar value will be ([Rauscher 



2012) 



f^new = ^exp + (^7th-t/cxp)(l-^,lab) . (11) 

Since X^iab depends on the plasma temperature, also the uncertainty factors are temperature de- 
pendent. 

For historical reasons, experimental uncertainties are usually given with linear error bars in- 
stead of uncertainty factors, although the latter would be more appropriate, especially for values 
approaching zero (suc h factors imply the use of logn ormally distributed errors instead of normal 



distributions, see, e.g.. iLimpert. Stahel. Abbtll200ll ). They can be mapped to uncertainty factors 



but a symmetric error ibA (with A > 0) of a value v (either a MACS or a reactivity) will lead to 
an asymmetric uncertainty 

V 

f^down =—!L^ . (12) 

In our definition uncertainty factors U > 1 and therefore the lower limit of v is given by v/U'^°'"'^ 
and the upper limit is given by vU^^. For symmetric uncertainty factors, as usually encountered 
in astrophysical rates, U = U^^ = [7'^°™" and v is in the range 

v/U<v<vU , (13) 

as used above. 

It is important to properly interpret the uncertainty factors defined above. Compatible defi- 
nitions of theoretical and experimental uncertainties have to be used especially in Equation (jlip . 
While experimental uncertainties are usually quoted as confidence intervals from stati stical dis- 



tribu tions, theoretical errors cannot be rigorously defined in such a manner (see, e.g., [Rauscher 



20121 . for a more detailed discussion). Often, theoretical upper and lower limits are given with the 
implicit assumption of a uniform probability for the actual value to be located anywhere within 
the limits. Moreover, experimental uncertainties contain statistical and systematical errors, which 
have to be treated differently. From the point of view of astrophysical application, i.e., variation of 
reaction rates, the range of values given by the uncertainty factors as shown in Equations (|lip and 
p3|) should cover a reasonable range of possible reaction rate values. In varying the rates, uniform 
probability is usually assumed within the given range, although Monte Carlo studie s also allow for 



normal and lognormal distributions (ISmith et al.l |2002| : ISmith &: NabereinevI |200J) . For practical 



purposes - although mathematically not rigorous - a simple but reasonable choice would be to 
assume a 2-sigma confidence interval for the experimental uncertainties Uexp and estimated upper 
and lower limits for theory values also containing 95% of the possible values, and then to assume 
that the actual values are uniformly distributed across both ranges. Using a 3-sigma interval of 
the (non-uniformly distributed) experimental data would be, in my opinion, exaggerated because 
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this would include outliers with low probability which are then assigned too high importance in the 
conversion to a uniform distribution. 

Since Xi is a theoretical quantity, it also has an inherent error. The impact of its uncertainty, 
however, is small with res pect to the expe r iment al and theoretical uncertainties in cross sections 



and rates. It was shown in lRauscher et al.l (120111 ) that the mamitude of the error scales inversely 



proportionally with the value ot X, i.e., X~has zero errol and that the nncertamty factor 

Ux > 1 of is given by max(nx, ^/ux), where 

ux = u{l-X,) + X, (14) 

and u is an averaged uncertainty factor in the predicted ratios of the Ri. In any case, the uncer- 
tainties are sufficiently small to preserve the magnitude of Xi, i.e., small Xi remain small within 
errors and large Xi remain large. The uncertainty factors Ux can easily be inserted in Equations 
(jlOP and (jlip and the uncertainties quoted for the reactivities in §[5] and Tabled] also contain Ux- 



5. Application to neutron capture reactions in KADoNiS 



The iKADoNiS v0.3l (|2009l ) compilation (see also iDillmann et al.l 120101 ) gives recommended 
MACS for neutron capture at 30 keV from experiment and their uncertainties. It also provides 
(n,7) reactivities in the range 5 < kT < 100 keV, derived either from experiment or from a 
combination of experiment and theory. Those reactivities did not consider the influence of Ajiab 
and therefore have to be revised. At the same time, updated error estimates can be provided. 

For cases of rates experimentally provided over the full temperature range, it is obvious that 
these are laboratory reactivities i?jiab, as discussed above. They are equal to the stellar reactivities 
R* as long as X^iab = 1. When X^iab < 1, a different approach has to be taken. For each temperature 
in the 5 < kT < 100 keV range, the theoretical stellar reactivity has to be renormalized, using 
Equation (llOp . Each value will then have an uncertainty according to Equation (jlip . 

The cases treated by a combination of experiment and theory were derived by renormalizing a 
theory value to an experimental value at a single temperature. The values at all other temperatures 
were then taken from the renormalized theory. A similar approach can be kept when including the 
influence of Ajiab . The theoretical stellar rate has to be renormalized according to Equation (|10p 
at one temperature, then values at all other temperatures are given by the renormalized stellar 
rate. Since the experimental uncertainty is only given at one temperature. Equation (jlip can be 
straightforwardly applied only when assuming that the experimental uncertainty factors would be 
the same at other energie^. 



^as long as G is known; this is the case close to stability 

very conservative approach - not adopted here - would include an experimental uncertainty factor only in the 
value for the measured energy /temperature and use pure theory uncertainties for all others. 
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What are the uncertainty factors C/gxp and JTth to be used in Equ ation (1110? The e xperi mental 
uncertainty is given for the recommended MACS at kT = 30 keV in iKADoNiS vO.a (|2009l ). It is 
stated in the compilation that the relative errors are expected to be similar at all other temperatures. 
As it is common in experimental nuclear physics, the quoted errors are on the 1-sigma level. 
Following the procedure recommended in § HI a 2-sigma error is used in the following to define an 
uncertainty range cov ering most of the allowed values. It is derived by doubling the error given in 
KADoNiS v0.3l ([20091) , and converted to C/exp by applying Equation (fT2]) . Since systematical errors 
are not given separately, they cannot be treated differently. 

The inherent theoretical uncertainty is difficult to estimate, as explained in ^ For neutron 
ca pture at stability the re i s a comparison to the compilati o n of iBao et al.l (j2000l ) (a predecessor 
of lKADoNiS vO.d 12009! ) in iRauscher. Thielemann, k Krat3 (jl997|). An average deviation of 30% 
is found, this w ould translate into an uncertainty factor Uth = 1-3. Comparing the MACS of 



KADoNiS v"a^ (12009.) to the theoretical (g.s.) MACS calculated with the code SMARAGD, vO.8.4 



(|E,auscheijl2009l . 



2011 



), an average deviation of 10% is found. On the other hand, local deviations 
up to a factor of 2.4 are present (but this may be due to a low level density, rendering the Hauser- 
Feshbach model inadequate). Moreover, the comparison is only made for MACS and not for stellar 
rates. If one wants to make sure that the uncertainty range of Equation (jl3p really covers the 
possible stellar values, it may be safer to assume C/th = 2. This is equivalent to a 2-sigma confidence 
interval, containing 95% of the possible values, were they statistically distributed. 

Finally, u has to be chosen to estimate the error in Xi. Most likely, 1 < u < Uth but variation 
of u within that range does not alter the results strongly. Therefore, u = C/th was chosen for 
simplicity. 

A re-evalu a tion of the stellar neutron capture rates for all cases with Z > W found in 
KADoNiS v0.3l (|2009l ) was performed using the above procedure. Th eoretical rates and MACS 



were calculated with the code SMARAGD, vO.8.4 (jRauscheijl2009l . l201ll ). The results are presented 



in Tabled] Each line initially specifies the element symbol, charge Z, and mass number A of the 
target nucleus. This is followed by 11 groups of values, for 11 temperatures kT: at 5, 10, 15, 20, 25, 
30, 40, 50, 60, 80, and 100 keV. Each group contains, in this order, X^iabCT), Ux{T), Kew(^)> and 
the final uncertainty factor U*^^{T) of -R^ew the given temperature. These uncertainty factors 
C/*g^(T) also define the range within which a stellar rate should be varied in reaction sensitivity 
studies. 



As already pointed out in IRauscher et al.l (1201 ll ). (n,7) reactivities in the higher mass part 
of the s-process are not well constrained by experiment. Figure [T] shows the uncertainty factors 
f7*g^ of the stellar reactivities at kT = 30 keV. Again, large uncertainties are found in this region 
already at that energy, despite of the fact that experimental uncertainties are small. The theoretical 
uncertainties are retained to a large extent. This is the consequence of a small laboratory level 
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contribution Xjiab in target nuclei with high nuclear level density, such as deformed nucleio At 
higher temperature, the uncertainties approach those of the model for all nuclei because Xjiab 
becomes small. Even high precision measurements cannot improve the uncertainty unless they are 
able to investigate transitions on excited levels separately and use theory to construct the stellar 
rate from this information, applying Equations 

An additional complication when performing renormalizations was mentioned in § 13.21 the 
question of reaction model applicability. This was ignored in the results presented here but will 
be treated in more detail in the upcoming KADoNiS version. It is expected, however, that the 
uncertainty stemming from the model (in) applicability is anticorrelated with Xi, i.e., the statistical 
Hauser-Feshbach model should be applicable when Xi is small because this implies a high level 
density. 

Furthermore, it has to be pointed out that also values for ^'''^Lu and ^^'^Ta are included in 
Tabled) It is not clear whether all levels in these nuclei are fully thermalized at the lower s-process 
temperatures and therefore these values may have a larger inherent error. Thermal population of 
all states is a fundamental assumption in the calculation of the stellar rates. 



6. Summary 

A general formalism to include experimental reaction cross sections into calculations of stellar 
rates was developed, which also allows to assess the reduction of uncertainties in the stellar rates 
by experiments. As an important exarnple fo r the application of the procedure, stellar neutron 
capture reactivities from lKADoNiS v0.3l (|2009l ) were revised and the remaining uncertainties have 
been shown. Although the uncertainties in the stellar rates are close to the experimental values for 
a number of nuclei, many of the uncerta inties remain consider ably larger. This has important con- 
sequences for s -process models (see, e.g.. lArlandini et al.lll999l ) and the interpretation of meteoritic 



data (see, e.g.. iQin et al. 



2008 



Burkhardt et al 



2OI2I ) because it allows the rates of some reactions 



to vary within a larger range than previously assumed. 



The revised reactivities will be included in the upcoming new version of KADoNiS. 



This research was supported in part by the European Commission within the FP7 ENSAR/THEXO 
project and by the EuroGENESIS Collaborative Research Programme. 



^There is one exception to this, the reaction ^^^0s(n j7) , for which neut ron transitions to excited states in 



Os 



were measured separately recently ( Mosconi et al. 2010a b; Fuiii et al. 2O10l ). Because of this additional information, 
the uncertainty on the stellar rate is similar to the one of the measurements. The current version 0.3 of KADoNiS, 
however, does not include these results. 
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Fig. 1. — Uncertainty factors t/*^^ for stellar (n,7) rates at kT = 30 keV. (A color version of this 
figure is available in the online journal.) 
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Table 1. Revised (11,7) reactivities R^^w their uncertainty factors U*^^ for 11 plasma temperatures T; laboratory level 
contributions to the stellar rate X^iab and their un certainty factors Ux a re also given. Negative values indicate that no 

information is available in lKADoNiS v0.3l ()2009l ) at that temperature. 



kT = 5 keV kT = 10 keV ... kT = 100 keV 

Element Z A X Ux C/*,„ X Ux R*c^ U*,„ ... X Ux R^w 

(cm'^mol^-'^s"-'^) (cm'^mol~-'^s~-'^) ... (cm'^mol^-'^s"-'^) 



w 


74 


182 


1.0000 


1.0000 


3.7772 X 10^ 


1.0620 


0.9997 


1.0003 


3.8496 X 10^ 


1.0626 


0.2142 


1.7858 


6.8995 X lO'' 


1.8875 


w 


74 


183 


0.9997 


1.0003 


9.4536 X 10^ 


1.0624 


0.9754 


1.0246 


8.7088 X lO'^ 


1.1069 


0.1504 


1.8496 


8.9071 X 10^ 


1.9237 


w 


74 


184 


1.0000 


1.0000 


3.5910 X 10^ 


1.0469 


0.9999 


1.0001 


3.2889 X lO'^ 


1.0471 


0.2411 


1.7589 


4.4374 X 10^ 


1.8694 


w 


74 


185 


0.9977 


1.0023 


1.1637 X 10® 


1.2253 


0.9746 


1.0254 


1.0062 X 10® 


1.2603 


0.2338 


1.7662 


5.1391 X lO'' 


1.8970 


w 


74 


186 


1.0000 


1.0000 


3.8086 X 10^ 


1.0829 


1.0000 


1.0000 


3.4755 X lO'^ 


1.0879 


0.2852 


1.7148 


2.9677 X 10^ 


1.8475 


Re 


75 


185 


1.0000 


1.0000 


2.5450 X 10** 


1.0879 


1.0000 


1.0000 


2.3289 X 10® 


1.0879 


0.6478 


1.3522 


2.1110 X 10® 


1.5630 


Re 


75 


186 


1.0000 


1.0000 


1.8500 X 10® 


1.1613 


0.9945 


1.0055 


1.8000 X 10® 


1.1613 


0.1179 


1.8821 


2.8200 X 10® 


1.1613 
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Note. — Table [T] is published in its entirety in the electronic edition of the Astrophysical Journal Letters. A portion is shown here for guidance regarding its form and content. 



